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We investigate laser-driven vibronic transitions of a single two-level atomic ion in harmonic and 
hard wall traps. In the Lamb-Dicke regime, for tuned or detuned lasers with respect to the internal 
frequency of the ion, and weak or strong laser intensities, the vibronic transitions occur at well iso- 
lated "Rabi Resonances" , where the detuning-adapted Rabi frequency coincides with the transition 
frequency between vibrational modes. These vibronic resonances are characterized as avoided cross- 
ings of the dressed levels (eigenvalues of the full Hamiltonian) . Their peculiarities due to symmetry 
constraints and trapping potential are also examined. 
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I. INTRODUCTION 

Laser cooled trapped ions have been much studied, 
both theoretically and experimentally, because of metro- 
logical applications as frequency standards, high preci- 
sion spectroscopy, or the prospects of realizing quantum 
information processing y|. The coupling between in- 
ternal and vibrational degrees of freedom of the ion in 
harmonic traps due to laser light has in particular been 
examined as a way to implement the Jaynes-Cummings 
model or other simple model Hamiltonians for quan- 
tum gate design, motional cooling or state preparation 

0, am ■ This coupling has been predominantly consid- 
ered for weak fields, £l R << lot, where lot is the trap's 
lowest (angular) frequency (for direction x, we assume 
that the two other frequencies in y and z directions are so 
large that the ion motion is effectively one dimensional), 
and il R is the optically "on-resonance" Rabi frequency, 

1. e., for the zero-detuning case in which the laser fre- 
quency lul and the ion transition frequency loq between 
two internal levels coincide, A := wj,- loq = |2g. In the 
laser adapted interaction picture, after applying the usual 
(optical) rotating wave approximation and dropping the 
constant term of the zero point energy, the Hamiltonian 
reads for a classical laser field, 
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r\ := k^xo is the Lamb-Dicke (LD) parameter, with xq 
being the characteristic length of the oscillator xo — 

, and fci the laser wavenumber; also a z — |e)(e| — 



\g)(g\, and a and a) are the annihilation and creation 
operators for the vibrational states {|n)},rt = 0,1,2,.... 
We are assuming that the internal states {\g), |e)} of the 



ion have an infinite lifetime, thus neglecting spontaneous 
emission. 

At this point it is customary to separate the Hamilto- 
nian as [f| 



H = H B + V(n R ), 

Hb = Htrap + Ha , 



(4) 
(5) 



to carry out a second interaction picture with respect to 
Hb = H(Q R = 0), i.e., the free motion of the electronic 
and vibrational degrees of freedom of the ion or "bare 
part". If the rapidly oscillating terms are neglected in 
a "second" or vibrational rotating wave approximation 
and, in addition, the Lamb-Dicke regime is assumed, 
rj << 1, three possible resonance conditions arise be- 
tween the detuning and the vibrational frequency. They 
define motional sidebands and state subspaces with sim- 
plified, effective Hamiltonians: 

• A = 0: "carrier" sideband, a purely internal reso- 
nance between bare states \g, n) and |e, n), without 
vibrational excitation. 

• A = — lot', first red sideband for vibronic transi- 
tions between \g,n+ 1) and \e,n) represented by a 
Jaynes-Cummings model Hamiltonian. 

• A = lot- first blue sideband for vibronic transitions 
between \g,n) and |e,n+l) (anti- Jaynes-Cummings 
model) . 

These transitions are the basis of the Cirac-Zoller and re- 
lated controlled-NOT quantum gates 0, IE • By choos- 
ing larger detunings, one can obtain in principle further 
vibronic resonances (motional sidebands) with non- linear 
couplings depending on higher orders of the Lamb-Dicke 
parameter, but they are not easy to realize in the labora- 
tory since efficient motional cooling is usually performed 
in the Lamb-Dicke regime 

The vibrational rotating wave approximation breaks 
down for larger £l R . This is but a particular example of 
the failure of the "secular approximations" for a strong 
perturbation [j| • There are however good reasons for ex- 
ploring also stronger fields since one of the consequences 
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of weak fields is that the vibronic transitions involving 
one vibrational quantum occur slowly: the effective Rabi 
frequencies for them scale with t]£Ir. Looking for faster 
quantum gates, Jonathan, Plenio and Knight 9] consid- 
ered the possibility of using instead stronger fields for 
the optically resonant case, A = 0, finding a different 
vibronic resonance condition, : = lot, for transitions 
involving one motional quantum and laser adapted inter- 
nal states |±) = (l/\/2)(|<7) ± |e)). These transitions had 
been already noticed by Moya-Cessa et al. 01 m a study 
of super-revivals, by means of a unitary transformation 
leading to the transformed Hamiltonian 

HQr 

Hmc — Htrap H 2~° z 

+ h-(<r- +*+) + — {a-a^(<j-+a+), (6) 
i){g\. Going to an inter- 



where cr_ = |.g)(e| and <j + - 
action picture with respect to H trap + ^^-(T z , and ap- 
plying again a rotating wave approximation to get rid 
of rapidly oscillating terms gives the vibronic resonance 
condition \Qr\ = lot and suppresses any dependence on 
the detuning, but, mirroring the previously commented 
failure of the RWA for "large" £Ir, this scheme is clearly 
not valid for sufficiently large detuning A, and requires 
A << lot- As pointed out in in the weak field limit 
the last term ( "spin flip term" ) in Eq. © does not go to 
zero so it does not provide a good perturbative scheme. 
Note also that higher, nonlinear resonances, of the type 
Qr = ktOT, fc = ±2,±3, are not immediately obvious 
using a second rotating wave approximation although, in 
principle, could be identified with a more sophisticated 
treatment 0. 

Aniello, Porzio, and Solimeno 11] could go beyond 
these two partial approaches, valid for weak fields or weak 
detunings, and set a more general resonance condition 
for the principal vibronic resonance, namely for Q = lot, 
where 



(^ + A 2 ) 1/2 >0 



(7) 



is the detuning-adapted Rabi frequency [23 ■ Their treat- 
ment, however, required a somewhat complicated series 
of transformations and relied on finding a zeroth order 
Hamiltonian diagonal in the bare basis of the bare Hamil- 
tonian Hb = H(CIr = 0). We shall present here a sim- 
pler treatment in which the initial Hamiltonian is sep- 
arated into the following zeroth order and perturbation 
terms 

H = H SD + W(ri), (8) 

where the "semidressed" Hamiltonian is given by H (77 = 
0), 



Hsd = hLo T a [ a —o z 



(| e )( 5 | + | 5 )(e|), (9) 



and represents the ion in the trap with a laser field which 
does not induce any motional coupling, whereas 

TiQr 



W{rf) 



(e"' (a+at) - l) \e)(g\+h.c 



(10) 



represents the coupling term. The explicit dependence in 
the argument of W(rj) emphasizes that H = Hsd when 
i] = 0. Notice that the leading term in Eq. QlOjl depends 
on t/Qr, whereas Hsd depends parametrically on lot, A 
and £Ir. This will determine the parameter domains ap- 
propriate for a perturbative scheme based on Eq. JHJ). 
All vibronic resonances discussed so far for weak field, or 
weak detuning will be identified directly from the degen- 
eracies of the "semidressed" eigenenergies of Hsd- With 
the coupling term W{rj), they become avoided crossings 
and their width and energy splitting are well character- 
ized by perturbation theory. In this way, the vibronic res- 
onance condition of Aniello et al. may be obtained and 
generalized from the crossings of the energy levels of a ze- 
roth order Hamiltonian in the whole range of £Ir and A 
values. Since the vibronic resonances correspond to Rabi 
frequencies equal to the vibrational transition frequen- 
cies, the term "Rabi resonance" , which we proposed pre- 
viously for a related waveguide confinement |l3| . seems 
appropriate for the general case. 

For small Lamb-Dicke parameters (the smallness cri- 
terion depends on the particular resonance as we shall 
see), the resonances are well isolated, which facilitates 
state manipulation and control. Other factors affecting 
any possible use of the resonances for information pro- 
cessing or state preparation are the speed of the transi- 
tion, determined by the effective Rabi frequency at the 
vibronic resonance, and its width, which will tell us the 
minimum stability required in the laser detuning and/or 
intensity. We shall see that both factors as well as the 
isolation of the resonances can be deduced from the en- 
ergy splitting, for which we shall provide explicit and 
resonance-dependent expressions. 

The formalism is also easy to apply for a hard-wall 
trap, and comparisons are made with the harmonic one. 
The hard-wall trap has been experimentally realized op- 
tically d4] or in a microelectronic chip version |l5j . and 
has generated recently much theoretical work in the field 
of ultracold atoms in low dimensions 0, 0, 0, 0, |2(J . 

In the next section we shall first rephrase the weak 
field case in a time independent framework following |j| . 
This will be a useful reference since the more general 
case of arbitrary laser intensity in Sec. Ill can be treated 
very similarly. In Section IV we compare the results for 
harmonic and square traps and the paper ends with a 
discussion and technical appendices. 



II. LOW INTENSITY FIELDS; TIME 
INDEPENDENT APPROACH 

In this section we present a time independent approach 
[3j which is essentially equivalent to the usual time- 
dependent treatment based on a second application of 
the RWA [JUl. For Q R < lo t HEl, it; is natural to re- 
gard the coupling between the laser and the ion as a small 
perturbation of the bare system consisting on the trapped 
ion, and decompose the Hamiltonian Q as in Eq. 
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A=-2o* T A=-t!U. A=0 A=to T A=2d) T 



FIG. 1: Bare (f2 = 0, dashed line) and dressed (Q/u>t = 0.2, 
solid line) energy levels in the resolved sideband regime as a 
function of the laser detuning. The Lamb-Dicke parameter is 
r\ — 0.4 in both cases. Note the degeneracies in the bare case 
at A = ±kujr become avoided crossings when the states are 
dressed (k th motional sidebands). 



ground and excited levels are connected by the perturba- 
tion and therefore do not cross. Moreover, these avoided 
crossings will be well-localized and isolated, because of 
the assumption of a low intensity laser. Good approx- 
imations to the dressed states at avoided crossings will 
be obtained by diagonalizing the 2x2 reduced matrix 
in the degenerate subspace. If we consider the crossing 
of the \g, n) and |e, n') levels at the A = En ' h E " = ku>T 
(k = n' — n) resonance, the matrix in the degenerate 
subspace spanned by \g, n) and |e,n') reads 

*Wi.,«o> = \ ( nl n T ) • ( 17 ) 

This operator governs the dynamics of the system at the 
resonance, leading to Rabi Oscillations with a frequency 
|f2„„'| between bare states \g,n) <-> \e,n'). These reso- 
nances are nothing but the k th motional sidebands men- 
tioned before . 

The eigenvalues of the matrix H17|) provide the energy 
splitting in the avoided crossing, 



H = Hb + V(Q R ). The bare part, H B = H(Q R = 0), is 
diagonal in the bare basis {\g, n), |e, n)} with eigenvalues 



hA 
hA 



(ii) 

(12) 



where E n — nhcor are the energy levels of the oscilla- 
tor. These bare energy levels cross each other whenever 
the detuning is on-resonance with the vibrational-level 
spacing, see Fig. ^ (dashed line), 



A = kuiT, 



(13) 



with k — 0, ±1,±2, .... When the laser is applied, 
these degeneracies are removed and the crossings become 
avoided crossings, see Fig. ^ (solid line). This is just a 
particular case of the "no level crossing theorem" , which 
states that a pair of energy levels connected by a pertur- 
bation do not cross as the strength of the perturbation 
is varied. The matrix elements of the coupling potential 



(g,n\V(n R )\e,n') = ~m* 
(e,n\V(n R )\g,n') = hn„ 



(14) 
(15) 



(the rest of terms are zero) , where the effective Rabi fre- 
quency O n „/ , defined by 



Or. 



n R {n\e^ a+a ^\n') = il n > 1 



(16) 



is the coupling strength between the motional levels n 
and n', calculated in Appendix [S] These matrix ele- 
ments are in general different from zero, so that all the 



Ar 



= h\n r 



(18) 



In leading order (see Appendix 1X1) . Ae nn / cx n R ri^ n ~ n '. 
A criterion for resonance isolation is Ae nn i << Klot, 
since Hujt is the energy difference between consecutive 
resonances. This leads to the condition £l R << lot for 
the carrier, and f]Q. R << lot for the first sidebands, which 
leaves some room for a relatively high r\ beyond the LD 
regime 0- 



III. ARBITRARY INTENSITY FIELDS 

If more intense lasers are to be used, e.g. for £l R > lot, 
the resonances described before are not isolated and new 
transitions with different sidebands occur. It is possi- 
ble however, to find well isolated avoided crossings for 
higher intensity lasers if we work in the Lamb-Dicke 
regime where it is more natural to regard the LD pa- 
rameter r\ as the perturbative parameter, as in Eq. JSJ, 
H = H$d + W(r]), instead of the on-resonance Rabi fre- 
quency used in the decomposition (0J. 

The semi dressed Hamiltonian H$ d is easily diagonal- 
ized, with eigenvalues and corresponding energy eigen- 
statcs 



|e».±) = 



E n ± 



2 ' 

i /A±n 



(19) 

g,n) + \e,n) I , (20) 



where the N± are dimensionless normalization factors, 



(A±ny 



2Q 



(fi ± A) 



(21) 
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These semidressed energy levels are degenerate at "Rabi 
Resonances", where the detuning-adapted Rabi fre- 
quency is resonant with some vibrational transition, 



\n' - n\ujT, 



(22) 



or Q% + A 2 



(n' — n) 2 Lo\. This also provides the 
motional sideband resonances in the low intensity limit 
<C lot, A): the semidressed states in Eq. (l2Tfl) reduce 
to the bare states and are resonant at A = (n' — ji)lot- A 
peculiarity of the present approach compared to the one 
in Sec. II is that, since SI is defined as a positive number, 
red and blue sidebands are formally treated in the same 
manner, see also ^lj pflj . 

At Rabi Resonances, the semidressed states involved in 
the crossing are coupled. The strength of this coupling 
is given by 



n,s\W(r,)\e. 
1 



(A + s'fl)(n\e 



irj(a-\-a 



2 VNlNt 
+ (A + sf7)(n|e-"'( Q+at > - l|n') 

hA 



if s 



l\n>) 



(23) 




if s ^ s' 



where s — ±1 is a shorthand notation to represent the 
"sign" of each state and 



C n 

Sn 



(n\ cosr](a + a')\n'}, 
(n\ sin 77(0 + cr )|n'). 



(24) 
(25) 



From parity arguments it is clear that, generically, all 
these coupling elements may be different from zero, con- 
necting all the semidressed states among them. For the 
A = case, however, only states with different vibra- 
tional parity will be connected, see Fig. where only 
alternate crossings are avoided. 

The eigenvalues of the 2x2 reduced matrix in the 
degenerate subspace, 



W n 



(€ n , + \W\e n ,+) 
(e n ,JW\e n>+ ) 



(e n .+ \W\e n , t - 
(e„,,_|W|e„ v 



(26) 



provide the energy splitting in the avoided crossing, 



Ae r , 



2n 



(27) 



x yj n 2 R (cw -c nn ) 2 +4 (A^cl n ,+n-'s 

with leading order in 77 (see Appendix 0) 



2 

nti ' 



A €ji fi 



+ 1 



HQ rij 



ly./ni 



— if / even 



1 if I odd 




^=4co T 



FIG. 2: Semidressed (j) = 0, dashed line) and dressed (77 = 
0.1, solid line) energy levels for the resonant A = case as a 
function of the Rabi frequency Q. Note the degeneracies in 
the semidressed case at Q — \n' — u\ujt are only removed when 
dressed for odd values of k = n' — n. Degeneracies surrounded 
by squares are real crossings of the energy levels, while the 
circles represent avoided crossings (which are not evident in 
the scale of the figure). 



In particular, Ae n! „+i = UrjClji y/n+ 1. Keep in mind 
that Qr in these expressions has to satisfy the resonance 
condition (J22J- 

The criterion for isolated resonances is Ae„ jn +; << 
Tiujt whereas the resonance width with respect to fl can 
be estimated from the slopes in Eq. (|19|l to be given also 
by the energy splitting. 

All vibronic transitions in Sec. II can be described in 
this manner, in particular we recover the splitting of Eq. 
(|18|l for A >> flu in leading order. At the sideband- 
motion resonances for k = ±1,±2, A is a multiple of 
lot so this condition is well satisfied. The energy split- 
ting of the carrier sideband, a purely internal transition, 
is also recovered in that limit, but it is clear that this 
resonance does not require A >> £Ir. In fact in Section 
II, the expression for the carrier splitting is obtained for 
A = and a small, finite £Ir << lot- Notice that in 
this regime H (77 = 0) already provides the carrier transi- 
tion splitting, so this regime lies outside the range of our 
perturbative scheme based on W(rj) breaking the degen- 
eracy. 

For the first vibronic resonance n = 0,n' = latA = 
the resonance condition is = lot which implies a 

splitting Htjlot and an isolation condition 77 << 1 equiv- 
alent to the LD condition. The isolation condition is 
therefore more stringent on 77 for these resonances than 
for the motional weak-field sidebands described in Sec. 
II. 

Another difference is the dependence of the effective 
Rabi frequency and the isolation condition on lot- Let's 



(28) recall that, for a given laser wavenumber, 77 cx lo t 



1/2 



Thus flRTj << lot is satisfied better and better for in- 
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FIG. 3: Logarithm of the ratio R := \\ 1 1+( /f2 00+( | for the 
first few sidebands (starting from the ground state) and for 
different values of the LD parameter r\. Note that for the 
first and second sidebands the coupling strength is larger in 
the harmonic potential, while for higher order multi-phonon 
transitions, coupling is stronger in the hard-wall trap. 



creasing lot, i-e., for tighter confinement irrespective of 
the type of resonance: for small Q R as in Sec. II, or for 
small A. However, the splitting, and thus the speed of the 
transition, decreases with increasing lot in the first case, 
whereas it increases in the second, because £Ir ~ lot- 



The crossing of the semidressed energy levels define the 
resonances as in the harmonically trapped ion. In the 
case of a hard-wall trap, the Rabi Resonance condition 
is given by 



\E n — E n i I 



hn 2 
2ma 2 



\n — n \ 



(31) 



The energy splitting in a given Rabi Resonance will 
be given by the same formal expressions as for the har- 
monic case, replacing Q nn i by Xnn'- In leading order in 
r?, defined now as k^ajix, 



ASl R n(n + I) 
l 2 (2n + l) 2 



„2A 



if I even 



2? if I odd 



(32) 



which should be compared with Eq. (|28|l . Different de- 
pendencies on the LD parameter r\ are observed. Note 
that for zeroth, first and second sidebands (I = 0, 1,2), 
both splittings are of the same order, while for higher 
order multi-phonon transitions (/ > 2), transitions in 
the rigid wall potential become faster if a strong con- 
finement is assumed (LD regime), see Fig. [3J In fact 
numerical comparison shows that the carrier transition 
is also stronger in a hard wall for typical values of the 
LD parameter, see Fig. 0] 



DISCUSSION 



IV. HARD- WALL TRAP 

We now consider the confinement of the ion in a hard- 
wall trap 



V(x) = 



if |x| < § 



if \x\ > § 



The treatment is essentially the same as for the harmonic 
trap, replacing exp[z?y(a + a^)] by exp(ifcix), and the en- 
ergies and eigenstates of the harmonic oscillator by those 
of the square potential. Now the motional energies 



E n = 



h 2 /n7T\ 2 



2m 



(29) 



with n = 1,2,... are no longer equally spaced. Note 
that the ground state in this trap is labelled with n = 1 
in contrast to the ground n = state of the harmonic 
potential. The hard-wall coupling strengths between the 
different motional levels (playing the role of £l n n' in the 
harmonic wall) are 

Xnn' = n R ((f) n \e lkLX \(f) n <) = fl R ; / 2 dx<f) n (x)<j) n ,{x)e ikLX , 



which are explicitly calculated in Appendix iBl 



(30) 



The perturbative schemes based on different decom- 
positions of the Hamiltonian H in Eq. (JIJ , to treat the 
effects of a laser on an ion in a harmonic trap, lead natu- 
rally to the identification of different families of vibronic 
resonances, i.e., transitions localized in a relatively small 
range of some varying parameter (such as the detuning, 
the laser intensity, a combination of both, or even the 
Lamb-Dicke parameter 25] ) in which the populations 
of both vibrational and internal levels oscillate. A well 
known family corresponds to the "motional sidebands" 
which appear as avoided crossings in the dressed level 
structure with respect to laser-ion detuning for weak- 
fields. Less studied are the vibronic resonances for zero 
detuning and varying laser intensity 0, 0, 0|. They 
may also be identified as avoided crossings with respect 
to the optically on-resonance Rabi frequency [l3| . 

It is in fact possible to identify and characterize all 
the above resonances with the same perturbative scheme 
|ll| . including moreover generalized resonances with re- 
spect to the detuning-adapted Rabi frequency. Our first 
contribution in this paper has been to simplify the treat- 
ment originally proposed by Annicllo et al. |ll|. Our ap- 
proach is based on a straightforward decomposition of the 
Hamiltonian Q into a zeroth order, "semidressed" term 
corresponding to H with zero Lamb-Dicke parameter, 
and a coupling term responsible for vibrational coupling. 
The resonances are at the degeneracies of a semidressed 
Hamiltonian (in the {A, f2#}-plane on concentric circles 
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APPENDIX A: COUPLING STRENGTHS IN THE 
HARMONIC POTENTIAL 

The coupling strength, or effective Rabi frequency, is 
given by 

ft„n' = n n/n = n R (n\e^ a+a ^\n') 

= Q R / dx<p n (x)<p n ,{x)e lkLX , (Af) 



where the ip n (x) arc the normalized cigenfunctions of the 
harmonic potential. This integral gives |22^ 



FIG. 4: Relative coupling strenght for a carrier \g,n g ) 
|e,n s ) transition in a harmonic (n g = 0, solid line) and in a 
hard- wall trap (n g — 1, dashed line) as a function of the LD 
parameter r\. 



where the detuning-adapted Rabi frequency is a multiple 
of the vibrational quantum) , and the energy splittings at 
the avoided crossings of the dressed levels are estimated 
from degenerate perturbation theory. These splittings 
determine the resonance isolation, its width, and the os- 
cillation frequency that gives the transition speed, all be- 
ing crucial factors in applications for quantum gates or 
state preparation. 

The formalism is very easy to adapt to hard- wall traps 
so that the results for harmonic and square traps can be 
compared. This comparison is our second main contribu- 
tion. The main result is that, to leading order in r], the 
carrier transition, and vibronic resonances implying the 
interchange of three or more vibrational quanta are more 
effective (i.e., faster, more strongly coupled) for the hard 
wall trap. Explicit expressions have been provided and 
possible applications of this intriguing behaviour are left 
for future research. 



where n < (n > ) is the lesser (greater) of n and n' and L" 
are the generalized Laguerre functions, defined by |23| 



LS(X) =£)(-!) 



k I n + a \ X k 



k=0 



n-k I k\ 



The leading term is given by |24 



\n — n'\\ 



(A3) 



(A4) 



APPENDIX B: COUPLING STRENGTHS IN THE 
SQUARE POTENTIAL 



if \x\ < f 



oo if |x| > f 



For a hard wall trap of width a, 
V(x)-- 

the coupling strengths are given by the integral 

Xnn' = Sl R {<t>n\e- ikLX \(l>n') = ^r[] dx(/> n (x)<j) n , (x)e ikLX , 

(Bl) 

where 4> n (x) are the normalized eigenfunctions of the in- 
finite well potential, 
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<f) n {x) = (x\(j) n ) = 



cos —x if n odd 

a a 



sin —x if n even 

a a 



satisfying the boundary conditions 4> n (— f ) = 0n(f ) = 0. 
This integral gives 



i_2 | sin if n — n' even 



D 



cos ^ if n — n' odd 



where D is a common denominator given by 
D = [(k L a) 2 - (n' - n) 2 ir 2 ] [(k L a) 2 - (n' + n) 
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Note that, as in the harmonic case, these coupling 
strenghts are real when connecting motional levels of the 
same parity, and purely imaginary otherwise. 

We define the LD parameter for this type of poten- 
tial in analogy with the harmonic case, in which the LD 
parameter could be written as 77 = with ph being 

the quasimomentum in Tvujt = |r5r- For a hard- wall 
potential 77 can be written as r> = with ph w being 

Phw 

the quasi-momentum unit in the energy eigenstates, i.e., 
Phw = Jin /a, so that the LD parameter is defined as 
77 := k^ajn. In terms of this LD parameter the coupling 
strength reads 



with 



D=[ V 2 - {n 1 - n) 2 ] [r, 2 - (n' + n) 2 ] . 

In the LD regime, the coupling strength to leading or- 
der in 77 is given by 



Xnn' ~ 



Ann'v 



(71 — n') 2 (n + n') 2 



rj 2 if n — n' even 
2t7/7t if n — n' odd 



= -r]Q R i 



\n~n'\'_ 



if n 



n even 



ttD 



cos ^ if n — n' odd 
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